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Abstract 

An analytical/numerical method has been developed to predict the static thrust 
performance of non-axisymmetric two-dimensional convergent-divergent exhaust nozzles. 
Thermodynamic nozzle performance effects due to over- and underexpansion are modeled 
using one-dimensional compressible flow theory. Boundary layer development and skin 
friction losses are calculated using an approximate integral momentum method based on the 
classic Kdrmdn-Polhausen solution. Angularity effects are included with these two models 
in a computational Mozzle Performance Analysis Code, NPAC. In four different case 
studies, results from NPAC are compared to experimental data obtained from subscale 
nozzle testing to demonstrate the capabilities and limitations of the NPAC method. In . 
several cases, the NPAC prediction matched experimental gross thrust efficiency data to 
within 0.1 percent at a design NPR, and to within 0.5 percent at off-design conditions. 
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Introduction 

Supersonic cruise transport aircraft and modem fighter aircraft with supersonic cruise 
and dash capabilities use convergent-divergent (CD) exhaust nozzles to attain supersonic jet 
velocities necessary for flight in the high speed regime. Over the past several decades, 
propulsion nozzle research has led to the development of multi-mission supersonic exhaust 
nozzles that can provide efficient operation and meet varying requirements over an extended 
flight regime. For example, at supersonic cmise, where the ratio of lift to drag is low and 
fuel consumption is high [1], the exhaust nozzle must adjust its area ratio to fully expand 
engine flow and provide the optimum external geometry for low drag. For takeoff, climb 
to cruise, and air to air combat maneuvers, the exhaust nozzle will need to adjust for rapid 
changes in engine pressure ratio and accommodate high mass flow afterburning settings. 
Finally, the nozzle may be required to vector or reverse thrust, suppress jet noise, or 
maintain a low-observable, compact profile with a minimum effect on radar cross section 
and low IR emission. 

Because of these diverse requirements, exhaust nozzle design is an integrated, multi- 
faceted effort, and is one of the most important technologies involved in the development of 
a supersonic aircraft. As such, an understanding of the flow physics and performance 
characteristics of CD exhaust nozzles is critical from experimental, theoretical, and 
computational standpoints. To meet this need, ongoing research efforts at NASA Langley 
Research Center are directed at generating an experimental database for new nozzle 
concepts [2], and developing, validating, and applying computational and analytical 
methods for the prediction of nozzle performance. 

This paper provides a concise discussion of the basic flow physics, off-design 
operation, and thrust performance characteristics of convergent-divergent exhaust nozzles. 
As part of this discussion, a detailed thermodynamic thrust performance model is 
developed over the practical range of nozzle operation, and an approximate integral 
momentum boundary layer method is derived from first principles to calculate boundary- 
layer development and skin friction losses. These models were combined with geometric 
loss estimates to develop a computational thrust performance prediction method for two- 
dimensional CD nozzles. A comparison of computational results with e.xperimental data is 
presented for several test cases and gives valuable insights into the nozzle performance 
characteristics discussed as well as illustrating the applications and limits of the prediction 
method. 


1 



Nomenclature 


S ymbols 

a Speed of sound 

a Nozzle flap divergence angle 

A ID flow area 

A* ID sonic flow area 

A9/A8 Nozzle expansion ratio 

B Local nozzle wetted perimeter 

Cf Flat plate skin friction coefficient 

Cft] Thrust efficiency coefficient, CFri=F/Fi 

CpG Gross thrust efficiency coefficient 

Cp Specific heat at constant pressure 

Cv Specific heat at constant volume 

6 Boundary layer thickness 

6* Boundary layer displacement thickness 

D Internal nozzle skin friction drag 

ACfg Gross thrust efficiency loss decrement, ACfg^ 
f Boundary layer development function, f=f(x) 

F Fixed geometry nozzle thrust 

Fi Fully expanded, ideal variable geometry nozzle thrust 
Fci Ideal convergent nozzle thrust 

F Normalized thrust, F=F/Po7Ag 

Nozzle divergence angular coordinate 
g Boundary layer development function, g=g(x, NPR) 

7 Specific heat ratio, y=C^Cy 

H] Boundary layer shape factor, Hi= 5*/0 

H2 Boundary layer shape factor, H2=5/0 

X Boundary layer development function, X^X<x) 

m Mass flow rate 

p. Dynamic viscosity at static temperature 

po Dynamic viscosity at stagnation temperature 
M Mach number 

M Boundary layer Mach number 

M* Sonic Mach number, M*=l 

V Kinematic viscosity 


2 



P Frcestream static pressure 

Pq Stagnation pressure 

0 Boundary layer momentum thickness 

© Power function of 0, 0=05/4 

r Boundary layer temperature recovery factor 

R Ideal gas constant, R=Cp-Cv 

Re Reynolds number 

p Density 

p Boundary layer density 

p* Sonic density 

AS Entropy rise 

a Skin friction coefficient compressibility factor 

To Flat plate shear stress 

T Freestream static temperature, absolute scale 

To Stagnation temperature, absolute scale 

T Boundary layer temperature, absolute scale 

u Boundary layer inner velocity 

U ID flow velocity, boundary layer outer velocity 

U* Sonic velocity 

CO Boundary layer viscosity power law exponent 

X Curvilinear wall coordinate 

y Wall normal curvilinear coordinate 

Subscripts 

1C First critical (Choke) 

2C Second critical 

D Design 

SE Shock at exit 

7, 8, 9 Nozzle reservoir, throat, exit stations 
oo Ambient 

a Angularity effect 

f Friction effect 

p Pressure 

V Momentum 

w Wall (y=0) 
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Abbreviations 


Axi Axisymmetric 

CD Convergent-Divergent 

NPAC Nozzle Performance Analysis Code 
NPR Nozzle Pressure Ratio, NPR=Po 7 /P«> 
ID One-Dimensional 

2D Two-Dimensional 
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Convergent-Divergent Nozzle Thrust Performance Characteristics 

The control volume used to calculate nozzle thrust is shown in Figure 1 . Written in the 
streamwise direction, the steady flow control volume equation for conservation of impulse 
and momentum states that 

X Forces = rh(Uout “ Ui„ ) ( 1 ) 

Figure 1: Nozzle Thrust Control Volume 

F 



For the control volume shown, with pressures written relative to ambient and zero 
freestream velocity (static conditions), 

F - (P9 - Poo)A 9 = mU9 

So, net thrust is a combination of momentum and pressure thrust at the nozzle exit: 

F = riiU9 + (P9 — Poo)A9 
= F +F 

Now, with mass flow at the nozzle exit defined as 

1119 = P9A9U9 

and with the introduction of the ideal gas law, 

P9 = P9RT9 

the nozzle mass flow is given by 

r p ^ 

*9 = ^9U9 


( 2 ) 

( 3 ) 

( 4 ) 

( 5 ) 

( 6 ) 
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Substitution of (6) into (3) nets 


F = 


^ P ^ 

^RT9y 


A9U9 +(P9 -Poo)A9 


But, 


U9 = M939 = 


and the thrust equation becomes 


Next, 


F = tP9A9M^+(P9-P„„)A9 


Pq = 


K^OIJ 


P„ = 


^P ^ 

oc 




07 


07 


(7) 


( 8 ) 


(9) 


( 10 ) 


( 11 ) 


so that 


Ao = 


'8 7 


F = (Pq7A8)- 


y D_ ^ 


^9 


v^07; 


VFo 7 Po7y 


When divided by Po 7 Ag, equation (13) yields a normalized thrust 


F = 


Pn7A 


07-^8 


Ag 


TM§ 




vFo7; 


P 9 P°c 
VFo7 Po7/ 


( 12 ) 


(13) 


(14) 


Now, with Station 7 at the reservoir, the following can be written 

P<» 1 


07 


NPR 


(15) 


where NPR is the nozzle pressure ratio, or the ratio of reservoir stagnation pressure to exit 
back pressure. Furthermore, 


P„ P 


07 


P07 P°o 


9 _ 


1 


NPR 


yPoo y 


(16) 
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With the substitution of these definitions and algebraic relations into (14), the control 
volume equation for normalized nozzle thrust looks as follows: 


F = 


^07^8 


^9 


NPR 


J 


P9 


NPR I 


-1 


(17) 


From this equation, the normalized momentum and pressure thrasts can be identified as: 


Fv _ Ag 

'yM? 

f"p "\1 

P 9 

Po7Ag Ag 

NPR 

P 


(18) 



Po7A 


07^8 


Ag 


I fP9 
npr|^p« 


\ 

-1 

y 


(19) 


Equations (17) - (19) can be written as functions of NPR and nozzle geometry alone, 
but for now these equations are in a good form to discuss the thrust performance 
characteristics of convergent-divergent (CD) nozzles. Figure 2 shows a plot of normalized 
nozzle thrust vs. nozzle pressure ratio for a typical fixed geometry CD nozzle and a variable 
geometry CD nozzle at static conditions. The variable geometry nozzle changes expansion 
ratio with NPR such that exit flow is always fully expanded (i.e., P 9 =Poo) and thrust is 
entirely due to momentum. The fixed geometry nozzle is shown to have a lower 
normalized static thmst over the entire NPR range except at the design point 
(NPR=NPRd). where P9=Poo and fully expansion is attained. 


Figure 2: Normalized Static Thrust vs. NPR 
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Figure 3: Nozzle Thrust Efficiency vs. NPR 



A measure of nozzle "thrust efficiency" can be obtained by dividing the normalized fixed 
geometry nozzle thrust ( F) by the normalized variable geometry nozzle thmst ( Fj ). Such a 

measure is shown in Figure 3, and more clearly illustrates the performance regimes of a 
fixed geometry CD nozzle. At the design nozzle pressure ratio, nozzle exit pressure is 
equal to ambient pressure, and flow is supersonic and fully expanded. At off-design 
conditions above the design pressure ratio, the CD nozzle enters the underexpanded regime 
where exit pressure is higher than ambient pressure. In this regime, the fixed geometry 
expansion ratio (A9/A8) is lower than the ideal expansion ratio. The nozzle is physically 
too small for full internal expansion, and flow must expand to ambient pressure outside the 
nozzle as shown in Figure 4. This external, unducted expansion corresponds to a loss of 
thrust efficiency. 

Figure 4: Underexpanded Nozzle Flow 
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Below the design pressure ratio, the CD nozzle enters the overexpanded regime. In this 
off-design regime, the fixed geometry expansion ratio is too large for complete expansion 
to occur at any given NPR. To satisfy pressure boundary conditions in the reservoir and at 
the nozzle exit, flow must recompress. At nozzle pressure ratios immediately below 
design, nozzle exit pressure is lower than ambient pressure, and recompression occurs 
outside the nozzle through a series of shock and expansion fan "cells", shown in Figure 5. 
This irreversible compression mechanism is responsible for the loss in thrust efficiency 
shown in Figure 3, and defines the off-design "externally overexpanded" regime. 


Figure 5: External Overexpansion 





Farther below the design pressure ratio, the nozzle expansion ratio is too large for even 
internal expansion to persist, and recompression occurs in the nozzle divergent section 
through a standing normal shock as shown in Figure 6. While this regime still contains 
overexpanded flow, overexpansion occurs inside the nozzle only; flow downstream of the 
normal shock is subsonic, and nozzle flow can adjust to ambient pressure at the nozzle exit. 
This regime is known as the internally overexpanded or "shock" regime, and occurs 
between the choke or first critical nozzle pressure ratio (NPRic) and the second critical 
nozzle pressure ratio (NPR2c)> where the shock is at the nozzle exit. 

Figure 6: Internal Overexpansion 
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As indicated in Figure 3, thrust efficiency decreases at a greater rate below the design 
NPR of a nozzle. For externally overexpanded flow (P 9 <Poo), equation (19) shows that 
pressure thrust is negative; i.e., low pressure acting on the nozzle exit creates a drag. In 
addition, pressures along downstream portions of the nozzle divergent walls will be less 
than ambient, creating an internal pressure drag. For the case of internally overexpanded 
shock regime flow, exit pressure equals ambient pressure, but pressures on the nozzle 
divergent walls are still lower than ambient and there is still an internal drag. More 
importantly, though, the nozzle also has a subsonic exit velocity and momentum thrust is 
significantly reduced. 

Figure 7: Critical NPR's vs. Expansion Ratio 



A9/A8 


Figure 7 shows the first and second critical pressure ratios and the design pressure ratio 
plotted against nozzle expansion ratio for the ID, adiabatic, inviscid flow of air with a 
specific heat ratio 7^1.4. The ID approximation shows that the shock regime occurs over 
a very small NPR range, and implies that the majority of thrust efficiency drop-off below 
the design NPR occurs in the externally overexpanded regime. In reality this is not usually 
the case, and the shock regime encompasses a significant part of the NPR range below 
design. 
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Growth of the shock regime is the result of shock induced boundary layer separation in 
the divergent section of the nozzle, illustrated in Figure 8. The predominant effect of this 
separation is a reduction in subsonic pressure recovery past the shock. Since exit pressure 
is fixed at ambient, the loss in pressure recovery results in higher than normal pressures 
inside the nozzle divergent section. The nozzle shock will compensate for this by moving 
upstream of its normal (unseparated) location and decreasing in strength. For every 
increasing NPR in the shock regime, then, separation causes a lag in nozzle performance 
by delaying downstream movement of the nozzle shock. This results in an extended shock 
regime and a shifted, higher second critical NPR. For convenience, the shifted second 
critical NPR will be denoted as the "shock at exit" NPR, or NPRse- R is important to note 
that the design NPR of a nozzle is independent of off-design, internal behavior, and so a 
fixed design NPR and an extended shock regime result in a smaller externally 
overexpanded regime. 

Figure 8: CD Nozzle with Shock Induced Boundary Layer Separation 



In most cases, the major effect of boundary layer separation in a fluid dynamic system is 
the resulting change in the system geometry caused by the separation. Translated to a CD 
nozzle, this implies that shock induced boundary layer separation delays operating points of 
a nozzle by changing the effective nozzle geometry. This is most evident in the extreme 
case of a nozzle with fully detached shock induced boundary layer separation, leading to an 
important analogy shown in Figure 9 which states that the separation point becomes the 
effective exit of a new, lower expansion ratio nozzle. At each NPR, then, the separated 
nozzle behaves like a shorter, lower expansion ratio nozzle with a shock at its exit. Thus, 
the nozzle will pass through the shock regime as a series of smaller nozzles operating at 
their "shock at exit" NPR's. The sum of the infinite number of individual smaller shifts 
forms a cumulative shift and determines the actual NPRse- 
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Figure 9: Separated CD Nozzle Expansion Ratio Analogy 




As discussed, overexpansion occurs because the nozzle expansion ratio is too large to 
maintain full expansion for a given NPR. In essence, then, separation becomes a "variable 
geometry" mechanism by which nozzle flow adjusts to a smaller, more correct expansion 
ratio. As might be expected, the most notable result of this separation geometry adjustment 
is an increase in static thrust efficiency, shown in Figure 10. This increase is the result of 
the natural tendency of an overexpanded fluid stream to reach a more efficient balance 
between internal energy and momentum, which translates into more efficient expansion and 
better conversion of total pressure and temperature into momentum and thrust. In addition, 
separation and the associated loss in pressure recovery past the separation point result in 
higher static pressures on the nozzle divergent flaps, reducing the internal drag discussed 
previously. In the limiting case of full detachment and zero pressure recovery past the 
shock, pressures all along the divergent walls will be equal to ambient pressure and the 
drag force will be eliminated. 

Figure 10; Separation Effects on Thrust Performance 


Thrust 

Efficiency 
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Though separation can be desirable for improving static thrust efficiency, the stability of 
shock - boundary layer interaction and corresponding separation is unpredictable in most 
cases, and this gives rise to serious design concerns. Asymmetric, unsteady, or transitory 
separation can have a detrimental or even catastrophic effect on almost any fluid machinery 
device, causing flow vectoring, pressure pulsations, and oscillatory flow patterns. In 
addition, though stable separation increases the static thrust performance of CD nozzles, 
wind-on aeropropulsive performance may suffer. Internal separation in a fixed geometry 
CD nozzle can increase the effective base area of an aircraft afterbody and couple with 
external separation to pump down aft facing surfaces, increasing afterbody drag. Thus, 
while static nozzle thrust efficiency increases, aeropropulsive "thrust minus drag" 
efficiency decreases, and a net performance decrease is realized. 


Gas Dynamic Modeling 

To write normalized thrust as a function of NPR, exit Mach number and exit pressure 
ratio (P9/P<x>) must be known as functions of NPR. To determine these variables, nozzle 
performance will be divided into two NPR ranges: 

Low End {NPRic<NPR^NPR2c} 

High End {NPR>NPR2c} 

Recall that the first critical pressure ratio (NPRic) is the point at which the nozzle chokes, 
and the second critical pressure ratio (NPR 2 c) is the limiting point in the internally 
overexpanded regime at which a shock is at the nozzle exit. This second critical pressure 
ratio refers to an "ideal" case, and not the separation shifted, "shock at exit" NPR. 


Proceeding with the analysis, the first and second critical pressure ratios can be 
determined with ID flow relations as follows. First, knowing P9=Poo at the first critical 
pressure ratio, a simple algebraic relation can be written 

^ = Z9_Zq9 /20) 

NPRic Pq7 Pq7 Po9 ^07 

Now, since nozzle flow is isentropic at the first critical pressure ratio, Po 9 =Po 7 , and 


1 P9 

NPRic Po9 


( 21 ) 
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In terms of local pressure ratio, the area - Mach number relation is [3] 


fv-lf 

2 ^ 

A _ 1 2 J 

U+u 


Y+i 

2(7-1) 


Y-l 


ll/2 


1 - 


i— 

Y 


P ^ 

vPoy 



^Pq> 


( 22 ) 


With AJA*=Ag/As, and after substitution of (21), the area - Mach number relation becomes 



Y-n 


1/2 f 




Y+l 

2(Y-I) 



y-l1 

1/2 


( ^ 1 

Y 


r 1 ] 

InpRicJ 



[npRicJ 


1 - 


(23) 


For air with a specific heat ratio y= 1.4, equation (23) can be solved for NPRic^l-89 to 
give the choke nozzle pressure ratio. Though not important at the moment, the design NPR 
of the nozzle (NPRd) can also be determined from (23), since the design point fulfills the 
isentropic requirements of the choke point math model (20, 21). For the design condition, 
however, (23) must be solved for NPRd> 1 89. This double solution corresponds to 
subsonic and supersonic solutions of the area - Mach number relation. 


Along these same lines, the second critical pressure ratio can be determined by coupling 
normal shock relations with the current math model. By denoting static and total pressure 
upstream of the shock as P and Pq respectively, and the same downstream conditions as P' 
and Pq', the following algebraic expression can be written 


P P P Pq 

Pq? P Pq Pq? 


(24) 


Now, flow up to the shock is isentropic, so Po=Po7- In addition, flow downstream of the 
shock is at ambient pressure, so P'=Poo. Keeping the static pressure ratio across the shock 
written as P'/P, and the static/total pressure ratio upstream of the shock as P/Pq, 


P 


07 


1 


r p Y p ^ 


Pq7 NPR2C 




( 1 ) 


(25) 
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The static/total pressure ratio upstream of the shock corresponds to supersonic flow at the 
nozzle exit for the second critical and all subsequent nozzle pressure ratios, since exit Mach 
number of the nozzle remains constant above NPR2C- The upshot of this relation is that 
the pressure ratio upstream of the shock is simply the inverse of the design NPR; i.e., 

^ = ^ (26) 

Po NPRd 

and equation (25) reduces to the following: 


NPR 2 C IPANPRd 


The static pressure ratio across the shock can be determined through a normal shock 
relation, where M is the Mach number upstream of the shock: 


P 7+1 VY+1 


The upstream Mach number can be determined easily using a ID relation [4] and 
equation (26). 


NPRd Po 


L 2 


M^ 


Solving (29) for M^ and substituting the resulting expression into (28) yields 


f _lzi 

^ 2 ^ { 1 \ y 

P "y + lU-lJtNPRD) 


-1 - 


Finally, substituting (30) into (27) and performing some algebra results in 


NPR2C = NPRd J — NPR J - 1 - 

(7 + l)(Y-l) D Y + 1 
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With the critical operating points of the nozzle defined, nozzle exit conditions can now 
be calculated, starting with the low end. The low end has already been defined as the 
internally overexpanded regime, and is characterized by a standing normal shock in the 
divergent section of the nozzle and subsonic velocity at the nozzle exit. In this regime, exit 
and ambient pressures are equal and thrust is entirely due to momentum. As such, the 
normalized static thrust equation reduces to the following; 


F = 


F 

^07^8 


A 9 

A^ NPR 


(32) 


So, for the low end, only the exit Mach number must be determined. Knowing the nozzle 
has a choked throat, it follows that 


Ag = A*7 


(33) 


so that 


(k V 




K^oij 


f K V 


vA^y 


vPo?y 


(34) 


For non-isentropic flow between any points I and 2 with no heat or work added [3] 


AS , 
— = In 
R 


^01 

vPo2y 


(35) 


or, 

(36) 
Pq2 

In addition, the Area-Mach Number Relation looks as follows [4]: 

Y+I 
2(Y-1) 

eAS/R 

(37) 


Ai M 2 


1 + ^^M2 
2 * 


By substituting (36) into (37), an adiabatic, non-isentropic, no work flow relation is 
obtained for area, Mach number, and total pressure between any two points in flow: 


A2 ^ Ml 
Aj M 2 


1 + 


1 - 1 - 


IzI 

Izl 

2 


Ml 

M? 


Y+l 

2(Y-I) 



(38) 
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Now, A 2 could be thought of as a fictitious sonic point of A 9 , and A| of A 7 ; i.e., 



< 

II 

< 

(39) 


> 

II 

> 

(40) 

From this definition. 

the following are true: 



M 2 = M 9 = 1 

(41) 


Ml = M 7 = 1 

(42) 

and. 

Poi =P07 

(43) 


P 02 = P09 

(44) 


When modified with these relations, equation (38) reduces to 


^9 _ Pq7 


or. 


A^ 


A7 = A9 


09 






(45) 


(46) 


So, by taking advantage of the sonic reference point concept, the area - Mach number 
relation (37) has been pared down to a simple expression relating sonic reference areas and 
total pressures. Now, (46) can be substituted directly into (34) which yields 


or, 


V 




vPo7 7 




07 


V D ^ 


V^9A^09A^07j 


f A. V 


v^8A*'o7; 




VA 97 


v»'o9; 


(47) 


(48) 


Noting that 


'9 _ 


07 


07 


NPR 


(49) 


Equation (48) becomes 


A9 

vAg; 


npr; 


p ^ 




(50) 
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For any point in compressible flow, isentropic ID relations give [3] 



( 51 ) 


Po 



Y-1 

2 





The product of the two relations in (5 1) and (52) looks as follows: 


Y+l 


f -1 

^ p ^ 

1 

2 

2(Y-1) 

1 + "'^ 

Ia* j 

vPq) 

” M 

Y + l 


L 2 J 


Substitution of (53) into (50) results in the following expression: 


Y+l 



f ‘ ' 

j_ 1 

2 

2 (Y- 1 ) 

1 + ’'^ ^ 


Inpr. 

' M9 

_Y + 1_ 


2 


which can be written as 


(52) 


(53) 


(54) 


1 

M 9 


1 +^ 

-1/2 

fA9l 

f 1] 

2 

2 



Inpr; 

Y + l 


Y+l 

2 (y- 1 ) 


Equation (55) can be rewritten as a quadratic equation for M<p- as follows: 


where 




A9 

f ^ 1 

2 

V Ag; 

InprJ 

Y + l 


Y+l 

2(y-]) 


(55) 


(56) 


(57) 


Using the quadratic formula, the solution of (56) for M 92 is 




1 

7-1 



( 58 ) 
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By performing some algebra and knowing that M 92 must be positive, the end result is 


M^=— !— 
Y -1 



2 (Y-1) 



(59) 


Equation (59) can be substituted into (32) to net the complete low end normalized thrust 
equation 


F = 


F 

^ 07^8 


A 9 

r Y ^ 

‘ 1 


Ag 

Iy-iJ 

'^NPRj 

V 


(60) 


This analysis was performed with no reference to shock location or strength; shock 
losses were written in terms of an entropy rise and represented by a total pressure - area - 
Mach number relation. Determining the shock strength and location is a simple matter from 
this point. Denoting total pressure upstream of the shock as Pq and downstream total 
pressure as Pq', the following can be written, since flow up to and past the shock is 
isentropic: 

P'o - Pq9 _ Pq9 P 9 _ Pq9 P°o _ Pq 9 r ^ ^ /gjN 

Po P07 P 9 P07 P 9 P07 P 9 ^NPrJ 

Substituting equation (52) for P 09 /P 9 yields 


Po 




1 

NPR 


] 


(62) 


Once the total pressure ratio across the shock is determined, using (59) for M 92 , the 
upstream Mach number can be found using a normal shock relation in the form [4] 


Po 


1 + ^^M^ 


-,JL 

Y-»r 




Y+1 Y+1 


Y-l 


(63) 


Finally, using (51) with A*=As, the area at which the shock occurs is found as 

Y+l 


M 


2 


Y + l 


2(y-I) 


(64) 


With knowledge of the nozzle geometry, the streamwise location of the shock can be 
determined from the shock area. 
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Now that the low end has been covered in detail, attention can turn to evaluating the 
high end. In this NPR range, the exit Mach number is constant and is the same as that at 
the second critical pressure ratio upstream of the exit shock. So, the remaining unknown in 
the normalized static thrust equation ( 17 ) is the exit pressure ratio, P9/P00. This ratio can be 
determined as a simple function of NPR as follows. First, an algebraic relation can be 
written in the form 


Z9.^Ii.P07 (65) 

Poo Pq 7 Poo 


Now, P07/P00 is the NPR, and P9/P07 is simply the pressure ratio corresponding to the exit 
Mach number. However, from ( 26 ), this ratio is known to be the inverse of the design 
NPR. Modified with this information, the relation is 



1 

NPRd 


NPR = 


NPR 

NPRd 


and the normalized high end thrust equation becomes 


A 

F = 


1 

> 

YM9 

^ 1 1 Y 

Pq 7^8 ^8 

NPRd \ 

NPRd NPR J 


( 66 ) 


( 67 ) 


Equation ( 29 ) can be used to give exit Mach number as a function of the design NPR. 
When substituted into ( 67 ), the resulting normalized high end thrust equation looks as 
follows: 
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From this equation, normalized momentum and pressure thrust can be identified as: 
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Note that the normalized momentum thrust is simply a function of the nozzle geometry, 
i.e., A 9 /A 8 and NPRd, and does not change with NPR. Normalized pressure thrust, 
however, varies with NPR, subtracting from the momentum thrust for NPR<NPRd 
( overexpansion) and adding to the momentum thrust for NPR>NPRd (underexpansion). 


The next step in this analysis is to determine the thrust of an ideal, fully expanded, 
variable expansion ratio nozzle such that a fixed geometry thrust efficiency can be defined, 
where the thrust efficiency coefficient is simply the ratio of normalized thmst to normalized 
ideal thrust: 
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Going back to (17), with P 9 /Poo=l for fully expanded flow, the normalized thrust equation 
looks as follows, where F is replaced by F, to denote the ideal, variable geometry nozzle 
thrust: 
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(72) 


In this case, both the expansion ratio A 9 /A 8 and exit Mach number are functions of NPR. 
For air with y= 1.4, equation (23) can be used to determine the design expansion ratio for 
any NPR above 1.89; for NPR's below 1.89, no "on-design" CD nozzle exists, and that 
range must be addressed separately. For now, however, the relation in (23) can be 
substituted into (72) to give 
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Equation (29) can again be used to determine exit Mach number. When (29) is substituted 
into (73), the variable geometry normalized thrust equation becomes 



Because no "on-design" CD nozzle exists below NPR=1.89, a convergent nozzle must 
be used to determine an ideal nozzle thrust. In this case, the control volume model is still 
valid and simplified, since a convergent nozzle has an expansion ratio A 9 /A 8 =l. For the 
convergent nozzle, the normalized thrust equation looks as follows, where Fci denotes the 
ideal convergent nozzle thrust: 
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For air with a specific heat ratio of 7^=1. 4, a convergent nozzle chokes at NPR=1.89. 
Below NPR=1.89 there is always subsonic flow and exit and ambient pressures are equal. 
In light of this, the simple ID relation in (29) can be used to determine exit Mach number. 
Modified with this relation, the normalized convergent nozzle thrust equation becomes 
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Now that the equations governing nozzle thrust performance over the low and high end 
regimes have been derived, it is important to note that, in most cases, the second critical 
NPR is so low that the low end regime need not be considered from a modeling standpoint. 
In addition, shock induced separation effects in a real nozzle dramatically alter and extend 
the low end regime into the NPR range where thermodynamics are governed by the high 
end model. Thus, for most nozzle applications, nozzle performance can be represented by 
high end thrust equations alone. Exceptions to this generalization are high expansion ratio 
nozzles with shallow divergence angles. Such nozzles have a wide, separation free low 
end regime and the second critical pressure ratio may fall into a range of practical interest. 
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Nozzle Thrust Performance: Loss Effects 

The analysis performed thus far has resulted in a model for nozzle thrust efficiency 
based on the thermodynamics and gas dynamics of energy conversion in a fluid stream. 
Inefficiencies predicted by this model are not due to loss effects, but rather, are the physical 
thermodynamic result of expanding flow in a fixed geometry duct over a range of non-ideal 
pressure ratios. Since this model was developed assuming the ID, adiabatic, isentropic 
flow of an ideal gas, real life departures from this ideal model are thereby due to non- ID 
and secondaiy flow effects, heat conduction effects, frictional and boundary layer effects, 
geometry induced loss mechanisms (internal wave radiation), and real gas effects. 

Over the entire range of nozzle operation, loss effects cannot be modeled in a simple 
fashion. However, at the design point of the nozzle, flow is fully expanded and most 
closely resembles that predicted by the ID model. At this condition, dominant loss effects 
can be readily and accurately predicted. As will be seen, losses calculated at the design 
NPR can be extended to other off-design, high end NPR's. Though such losses are not 
applicable to low end NPR's, this regime need not be modeled in most cases as discussed. 

To include loss effects in the thermodynamic thrust performance model, it is first 
necessary to break nozzle thrust efficiency down into its momentum and pressure portions. 


C 



p j-p p +p 

1 - = y. 

Fi Fi 


( 77 ) 


In a similar fashion, a gross thrust efficiency coefficient, that is, one that includes losses, 
can be defined, where the normalized momentum and pressure thrust with losses are 
denoted by primed quantities: 


CpG = 


F +F 

^ V ^ ^ p 


Fi 


( 78 ) 


At the fully expanded design point of the nozzle (i.e., NPR=NPRd) the following is true: 
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where the peak gross thrust efficiency, C^q, is the ideal peak thrust efficiency of 1.0 
minus losses at peak: 

CP^ = 1-ACP^ (82) 

With these relations, the thrust efficiency and gross thrust efficiency equations become: 
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Combined, these two relations give 
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that is, the momentum thrust with losses is simply the ideal momentum thrust multiplied by 
the peak gross thrust efficiency. This result is of great importance, for it means that losses 
need only be predicted at the design point of the nozzle to model performance in the high 
end. Physically, this implies that nozzle losses affect only the momentum portion of 
normalized thrust, which was shown previously to be independent of NPR in the high end. 
From the peak gross thrust efficiency then, pressure thrust effects determine the actual 
thrust efficiency of the nozzle at each high end NPR. Thus, the peak gross thrust 
efficiency becomes the effective maximum thrust efficiency of the nozzle from which the 
thermodynamic model can predict over- and underexpansion inefficiencies at off-design 
high end NPR's. With these results, the gross thmst efficiency coefficient can be re- 
written as: 
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For comparison between the thermodynamic model and experimental data, a measured 
peak gross thrust efficiency can be inserted into equation (86). For an a priori performance 
prediction, however, the peak gross thrust efficiency must be determined by estimating the 
appropriate loss effects. The simplest such loss effect is angularity, which results from the 
non- ID flow of fluid at the nozzle exit. Angularity losses occur when the nozzle geometry 
generates non-axial jet velocities, as in the case of a nozzle with a non-zero exit divergence 
angle. For a fixed total thrust, components of momentum thrust generated in non-axial 
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directions will reduce the axial thrust of the nozzle, and a decrease in thrust efficiency is 
realized. For a 2D nozzle with a divergence half angle of a, the effects of angularity can be 
determined as follows. 

Figure 11; Nozzle Divergence Polar Coordinate System 
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First, using a polar coordinate system originating from the vertex of the nozzle 
divergence angle as shown in Figure 1 1, the nozzle exit plane can be broken down into a 
series of streamtubes, each corresponding to a polar element d(l). Each streamtube can then 
be thought of as containing an element of thrust dFy of the total ideal momentum thrust Fy, 
where dFy is defined as follows: 
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The component of elemental momentum thrust in the axial direction, dFya, is simply 


dFya = dFy cos(j) = 
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It follows that the total axial thrust is simply the elemental axial thmsts integrated over the 
nozzle exit plane: 
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Now, the gross thrust efficiency decrement due to angularity is: 
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but at the design NPR, Fi=Fv , and the decrement at peak is: 

=1-^ = 1-— (2D Nozzle) (91) 

’ Fy a 

Though not necessary for the current analysis, a similar derivation for an axisymmetric 
nozzle gives: 

(Axi Nozzle) (92) 

The second loss effect that can be predicted in a nozzle is skin friction drag, though this 
loss effect requires a more rigorous analysis than angularity. The internal friction drag of 
the nozzle can be determined by integrating the product of the local skin friction coefficient, 
dynamic pressure, and wetted perimeter over the nozzle streamwise length run. 

L 

D= jBCf(ipU^)dx (93) 

0 

The loss in thrust efficiency due to frictional effects can be obtained by dividing this 
internal drag by the ideal nozzle thrust, Fi: 

ACFG,f=-^ (94) 

Including angularity and friction losses, the peak on-design gross thrust efficiency of 
the nozzle is predicted to be 

Cf^ = l-Acg;„-Acfef (95) 

Nozzle Boundary Layer Development Modeling 

In its current state of technology, the science of computational fluid dynamics has 
evolved to a level of maturity where CFD can be readily applied and results are widely 
accepted. The advent of modem super computers has resulted in the ability to obtain 
solutions to the full Navier Stokes Equations that were unimaginable just twenty years ago, 
leading to detailed boundary layer predictions and an almost unlimited source of flowfield 
information. While this capability greatly increases the depth of the researcher's toolbox, it 
also demands an increased commitment; CFD analysis can be time consuming in a first 
approach, and in many cases, justifies an separate research effort of its own. 
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For many fluid dynamic applications, methods based on the approximate solution of 
boundary layer integral conservation laws can yield valid results. These classic methods 
are simple in derivation, numerical implementation, and application, and have been shown 
to calculate boundary layer development and skin friction drag with a high degree of 
accuracy and reliability. In general, however, integral methods are limited to a class of 
problems in which flow physics are well understood and detailed flowfield information is 
not required. For basic nozzle geometries, the prediction of nozzle thrust performance is 
one problem which fits into this category. 

Integral methods are based on solutions of the boundary layer integral equation for 
conservation of momentum, derived physically by von Karmen and mathematically by 
Polhausen in 1921 . The simplest approximate solution to this equation is the well known 
Karm^-Polhausen Method, used to calculate laminar incompressible flows over a flat plate 
with zero pressure gradient. Though readily adapted to compressible cases with pressure 
gradient, laminar methods of this type are not applicable to nozzle flows where Reynolds 
Numbers are too high for laminar flow to persist, and intense turbulence results in 
transition Reynolds Numbers that are typically lower than the equivalent flat plate [5]. 

In 1951, Tucker [6] developed an integral approximation method for calculating the 
turbulent boundary layer development in compressible flows, and in 1954, Bartz [5] 
derived a method for calculating the turbulent boundary layer development in axisymmetric 
rocket nozzles. Both turbulent methods were based on an "analogy" with laminar 
boundary layer theory; that is, laminar conservation laws were used to derive the basic 
governing equations. Turbulent boundary layer relations (vJ\J=f(y/d), Cf, etc.) were then 
imposed on the laminar model to come up with an analogous turbulent boundary layer 
model. Though Tucker was able to verify results from this model with experimental 
boundary layer development data, neither author attempted to make or validate skin friction 
drag calculations, possibly due to a lack of suitable data. The following is an extension of 
the approximate integral momentum method to a non-axisymmetric, two-dimensional 
nozzle for the calculation of turbulent boundary layer development and internal performance 
losses due to skin friction drag. 
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The Karman-Polhausen integral boundary layer equation for steady flow conservation 
of x-momentum is presented in (96), where the x and y axes are the local directions along 
and normal to the solid boundary [7]. 



Note that in the following analysis, u is the boundary layer velocity, U is the freestream or 
edge velocity, p is the boundary layer density, p is the freestream density, and P is the 
freestream pressure, henceforth assumed constant in y through a laminar boundary layer 
analogy. As part of this model, flow outside the boundary layer is assumed isentropic and 
ID. For this flow, the Euler Equation for conservation of momentum along a streamline 
can be written as follows [8]: 
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The outer pressure gradient can be imposed on the boundary layer model by substituting 
(97) into the integral momentum equation: 



Combining terms under the integral results in 
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Next, using the product rule to make the following expansion. 


0pu _ 1 dpu^ ~ 3u 
dx u dx dx 

the bracketed terms in (99) become: 
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Now, use of the product rule to expand the first bracketed term in (101) results in 
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Substitution back into the integral momentum equation nets 
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Making use of the following definitions for the boundary layer momentum thickness, 
boundary layer displacement thickness, and wall shear stress [7]: 
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the boundary layer integral momentum equation reduces to: 
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C 

. ^ ipU2 


(107) 


(108) 


29 



the integral momentum equation can be written as: 
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Finally, since 0=0(x), U=U(x), and pU2=plj2(x), d/dx becomes d/dx and (109) is: 
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Equation (1 10) can be put into a simpler form by noting that, from ID gas dynamics, with 
M=U/a as the freestream Mach number [5], 
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Substitution results in 
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Results obtained up to this point are based on laminar boundary layer theory. To extend 
this model to an "analogous" turbulent boundary layer, a fully developed turbulent 
boundary layer with a 1/7-th power law similarity velocity profile (u/U=(y/5)l^) and a 
turbulent flat plate skin friction relationship can be assumed. The local nozzle skin friction 
coefficient can then be determined using a Blasius skin friction equation for a turbulent 
boundary layer. The Blasius equation given by Eckert [9] is: 
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For a compressible boundary layer, the dynamic viscosity in (114) is a function of 
temperature and must be calculated at the surface "wall" temperature. In the case of a heat 
conducting compressible boundary layer, this would require a simultaneous solution of the 
boundary layer integral energy equation. For cold nozzle flows, however, the simplifying 
and reasonable assumption of an adiabatic boundary layer can be made, since subscale 
nozzle testing is usually performed over a length of time such that temperature equilibrium 
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between flow and the nozzle wall can be attained. Noting that the freestream and boundary 
layer temperatures are T and T, respectively, and making the additional assumption of a 
similarity boundary layer temperature profile where T / T = /(u / U), the energy equation 

can be directly integrated to give [7] 

t = T + r^^[TM^- W] (115) 

where M = u / a is the boundary layer Mach number and r is the boundary layer 
temperature recovery factor, which accounts for temperature increases through viscous 
dissipation in the boundary layer. For a compressible turbulent boundary layer with a 
Prandtl number of 0.72, r is approximately 0.88 [7]. 

At y=0, M=u=0, T = f and (115) gives 
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To relate viscosity to temperature in a compressible flow, Schlicting [7] suggests the use of 
a simple power law: 
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where values of o>=0.75 showed good agreement with experimental data in the temperature 
range from 360 to 560 °R. From ID gas dynamics, local temperature is related to 
stagnation temperature by [3] 
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Substitution into the skin friction coefficient equation (114) gives 
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where 0 is defined as follows: 
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From the continuity equation. 


pUA = p*U*A^ 


A* 

A 


which gives the following: 

pU = p*U*l 

leading to the final form of the skin friction coefficient. 
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Substituting the compressible skin friction equation (124) into the integral momentum 
equation (113), multiplying that equation by and noting that 
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the result is: 
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Note that the boundary layer shape factors Hi=5*/0 and H2=5/0 are present in the final 
form of the integral momentum equation. These are evaluated as follows. First, using the 
ideal gas law with pressure constant through the boundary layer in y, the following can be 
written: 
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With this, equation (1 15) can be used to rewrite the momentum and displacement thickness 
integrals in (104) and (105) as: 
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and. 
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Finally, with u/U=(y/5)l/7, the integrals become: 
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Due to their complicated form, equations (130) and (131) are best solved using quadrature. 
For these equations, Simpson's Rule gives: 
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X X 

5* = -[c(0) + 4c(5 / 2) + c(6)] = - 
6 6 


With the following substitution 
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= 1 + 0. 1797 r| ( 1 34) 


the shape factors are: 
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Finally, with the following functional definitions, 

t. 

1-M 


f(x) = f 
4 


YM^-Hi-2 

r2 


Afta A 


dx V A* 


Vt' ^ 

FO 

1 ] 

f— V 


IH2J 



g(x,NPR) = -(0.0228)a| 


0 = 05/4 

the integral momentum equation can be written in a compact form. 
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Equation (140) is a first order, linear, nonhomogeneous ordinary differential equation 
which sets up an initial value problem when an onset boundary layer momentum thickness 
is specified. Several different numerical schemes have been used to solve this problem 
with a similar degree of success, including a multistep backward-space (upwind) finite 
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difference method, a fourth order Runga-Kutta method, and a fourth order Adams- 
Bashforth predictor - Adams-Moulton corrector method. 

As mentioned at the beginning of this analysis, the pressure gradient used to calculate 
the boundary layer development is based on the inviscid outer flow. Once an inviscid 
solution is obtained, the estimated boundary layer displacement thickness can be added to 
the inviscid wall geometry and the boundary layer calculation iterated upon until a solution 
is obtained that satisfies both the boundary layer model and the pressure gradients imposed 
by the boundary layer displacement. Pressure gradient changes due to displacement effects 
are generally on the order of 5* itself, so an iterated solution may or may not vary 
significantly from the initial, inviscid boundary solution, depending on the nozzle geometry 
and displacement thickness development. However, an iterative procedure of the type 
discussed is easy to implement and generally converges quite rapidly, so little additional 
computational work is required to correct for displacement effects. 

Code Development and Application 

The boundary layer integral momentum method derived above was combined with 
thrust performance calculations and integrated into a prediction program NPAC, which 
stands for Nozzle Performance Analysis £ode. Results presented in this paper were 
obtained using a 4th order Runga-Kutta method embedded in an iterative boundary layer 
displacement correction scheme to solve the integral momentum equation and calculate 
boundary layer development and skin friction drag. For this study, NPAC was run on a 
Digital Equipment Corporation Alpha 3000 workstation in interactive mode. 

The non-axisymmetric nozzle geometries used in this analysis and the data presented in 
this discussion originated from sub-scale nozzle testing conducted at the NASA Langley 
Research Center 16 Foot Transonic Tunnel Complex over the past 20 years. Propulsion 
test facilities at this complex use a high pressure air supply system with a series of 
reservoirs, valves, filters, and heat exchangers to provide clean, dry airflow at a constant 
total temperature of about 530 °R and mass flow rates up to 30 Ibm/sec at atmospheric back 
pressure. Details of a typical propulsion simulation system are shown in Figure 12. 

Nozzle air flow rates are calculated based on pressure and temperature measurements in a 
calibrated multiple critical venturi system located upstream of the propulsion simulation 
model. Forces and moments generated by the nozzle are measured by a six-component 
strain-gauge balance mounted on the propulsion model centerline and corrected to 
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calibration standards. Stagnation conditions are measured by pitot pressure rakes and 
thermocouple probes. Reference [10] provides further details of the test facilities at the 
16 Foot Transonic Tunnel Complex. 

Figure 12; Propulsion Simulation System 



For this discussion, four different nozzles were analyzed and results compared to 
experimental data to show both the validity and the limitations of NPAC. All nozzles were 
non-axisymmetric, two-dimensional, convergent-divergent nozzles with nominal throat 
areas (As) ranging from 2.50 to 4.32 in^ and expansion ratios (Ag/Ag) ranging from 1.09 
to 1.80. For analysis, the centerline length of each nozzle was broken down into 2000 
intervals of approximately 0.005 in. each in length. Each interval was projected to the 
nozzle channel contour to determine the local Runga Kutta step size (Axi) which was then 
halved to provide starting (xj-Axi), intermediate (xi-Ax,/2), and calculation (xi) points for 
the fourth order Runga Kutta solution. It is important to note that the Runga Kutta solution 
was largely unaffected by changes in interval size, and that size was chosen primarily to 
obtain proper resolution of the nozzle geometry for wetted area calculations. 

At the appropriate design NPR, boundary layer development calculations were started 
downstream of the propulsion simulation system choke plate with an assumed zero 
thickness onset boundary layer. Boundary layer displacement thickness corrections were 
made to the inviscid nozzle geometry until the boundary layer displacement thickness 
development converged to machine accuracy. At that point, skin friction drag calculations 
were made over the full length of the nozzle flaps and sidewalls beginning at the nozzle 
entrance. The skin friction estimate was combined with the geometric angularity loss 
calculation to predict the peak gross thrust efficiency of the test nozzle. From that point. 
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the thermodynamic thrust performance model was used to determine over- and 
underexpansion effects at off-design nozzle pressure ratios. 

Case 1 

The nozzle studied in Case 1 was tested by Mason, Putnam, and Re [1 1], and has a 
nominal throat area of 4.3106 in^, a divergence angle of 1.21°, an expansion ratio of 1.09, 
and a corresponding design NPR of 2.97. Geometry details are shown in Figure 13. This 
nozzle was used as part of an investigation to determine the effects of throat contouring on 
2D-CD nozzle static internal performance. 

Figure 13: Case 1 Nozzle Geometry. Dimensions in inches. 



The calculated boundary layer development for this nozzle is presented in Figure 14 
after one iteration (inviscid pressure gradient), and after the last iteration (displacement 
corrected pressure gradient). The only difference in the boundary layer development in the 
two cases is at the nozzle entrance, where the corrected pressure gradient resulted in a 
slightly thiimer boundary layer due to displacement effects in the subsonic "flat plate" 
instrumentation section upstream of the test nozzle. This initial difference is seen to be of 
no consequence further downstream, where the boundary layer development is the same 
for the inviscid and corrected cases. As observed by Bartz [5], the nozzle boundary layer 
reaches a minimum thickness slightly upstream of the nozzle throat and growth in the 
divergent section of the nozzle is nearly linear. Other nozzle cases investigated in this 
analysis displayed similar boundary layer development. 
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Figure 14: Calculated Turbulent Boundary Layer Development for Case 1. 



Iteration Number 
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Convergence of the boundary layer displacement thickness development under the 
iterative correction scheme is shown in Figure 15, where Loo and L 2 norms of the residual 
vector are plotted against iteration number. Note that the residual vector consists of 
components representing the residual at each calculation point along the nozzle length. 

The correction scheme is seen to converge to machine accuracy in less than 50 iterations. 
Other nozzle cases discussed in this paper displayed similar convergence behavior. 

The gross thrust efficiency prediction for the Case 1 nozzle is compared to experimental 
data in Figure 16 and shows excellent agreement throughout the NPR range tested. Peak 
thru.st efficiency was predicted to within 0. 1 percent, and the NPAC prediction agrees with 
all data to within the measurement system precision (about 0.5 percent). Losses predicted 
by NPAC at peak were as follows: ACfg,P= 0 01 1. ^CpG.a-O. 

Figure 16: Predicted and Experimental [11] Thrust Efficiency for Case 1 



NPR 


Case 2 

Experimental results for the 2D-CD nozzle investigated in Case 2 are discussed by 
Berrier and Re in [12]. As shown in Figure 17, this nozzle has a nominal throat area of 
4.3262 in2, a divergence angle of 5.38°, an expansion ratio of 1.25, and a corresponding 
design NPR of 4.22. This particular nozzle was used as part of a test program to evaluate 
the effects of several geometric parameters on nozzle internal performance. 
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Figure 17: Case 2 Nozzle Geometry. Dimensions in inches. 



For Case 2, the gross thrust efficiency prediction is compared to experimental data in 
Figure 18. Once again, there is good agreement throughout the NPR range shown, and 
peak thrust efficiency fell between repeat experimental data points at NPR=4.25. For this 
nozzle, NPAC predicted a peak gross thrust efficiency of 0.991, with a loss breakdown of 
ACpG^f=0.007 and ACpc cc=0.002. 

Figure 18: Predicted and Experimental [12] Thrust Efficiency for Case 2 
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Case 3 

The nozzle used in Case 3 has a nominal throat area of 4.3 172 in^, an expansion ratio of 

I. 797, a corresponding design pressure ratio of 8.78, and a steep divergent flap angle of 

I I. 01°. This nozzle was tested by Hunter [13] to investigate passive shock - boundary 
layer interaction control concepts. Specific details of the nozzle internal geometry are 
presented in Figure 19. 

Figure 19: Case 3 Nozzle Geometry. Dimensions in inches. 

U 2215 2.275 A 



Unlike previous cases, the thrust efficiency prediction for Case 3 in Figure 20 shows 
less favorable agreement with experimental data over the entire NPR range tested. NPAC 
did predict the peak gross thrust efficiency coefficient and modeled off-design performance 
almost exactly down to NPR= 5.0. Below this NPR, though, the measured gross thrust 
efficiency rose higher than the NPAC prediction, and the difference increased at lower 
NPR's. Such behavior might seem puzzling at first, since losses measured at the design 
NPR should not change at off-design high end NPR's, and would certainly not decrease at 
the rate seen in experimental data. As discussed earlier, however, one mechanism can 
significantly alter off-design nozzle performance and would increase thrust efficiency: 
boundary layer separation. Experimental off-design internal pressure data for the Case 3 
nozzle are presented in Figure 2 1 , and show strong evidence of shock induced separation 
with little or no downstream static pressure recovery in the divergent section of the nozzle 
for NPR's between 1.8 and 5.0. Above NPR=5.0, the nozzle was shock free, and thus 
was performing as the NPAC model would predict. 
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It is interesting to note that between NPR's 1.8 and 2.4, experimental gross thrust 
efficiency data followed the same trend as the NPAC prediction, though it was shifted to 
lower NPR's by ANPR~1.0. This supports the notion of an NPR range shift discussed 
previously, and indicates that a large portion of that shift occurs at lower NPR’s, where the 
nozzle shock is further upstream and more dramatic changes in the effective nozzle 
geometry can occur through separation. 

Note that for Case 3, the predicted gross thrust efficiency at peak was 0.986, with 
associated losses of ACFG,f=0-008 and ACFG,a=0-006. 

Case 4 

The final nozzle studied in this analysis was tested by Capone and Berrier [14] as part of 
a wind tunnel experiment on a 1/10-scaIe, twin engine F-18 prototype aircraft model, and is 
a 2D-CD nozzle with an expansion ratio of 1.15, a design NPR of 3.46, and a nominal 
throat area of 2.50 in^. Details of the nozzle geometry are presented in Figure 22. Note 
that this nozzle had a cutback outside sidewall and an extended nozzle inter-fairing on the 
inside sidewall, and thus was not truly a 2D nozzle along its entire length. The geometry 
shown in Figure 22 is defined by the nozzle upper and lower flaps and does not take 
sidewall variations into account. 

Figure 22: Case 4 Nozzle Geometry. Dimensions in inches. 
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The thrust efficiency prediction for this nozzle is shown compared to experimental data 
in Figure 23. In this last case, the NPAC prediction does not match experimental data very 
well. Most of the thrust efficiency data at and below design is within a 0.5 percent 
measurement system precision band around the NPAC prediction, but previous cases have 
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demonstrated that agreement should be much better than that quoted limit. In addition, 
experimental gross thrust efficiency was much higher than the NPAC prediction at 
underexpanded NPR's. Regardless of exact levers, However, there is a notable NPR shift 
between the NPAC thrust performance curve and experimental data over the entire NPR 
range tested, and a closer look at the data shows that the nozzle reached a peak thrust 
performance level at about NPR=4.5, approximately 1 NPR higher than the ideal peak. It 
is apparent that physics outside of the realm of the NPAC gas dynamics / integral 
momentum boundary layer model were present in this case. 

Figure 23: Predicted and Experimental [14] Thrust Efficiency for Case 4 



The flow physics missing from the NPAC model in this last case can be linked to two 
relatively simple factors: the cutback outside sidewall, and the twin engine configuration. 
The cutback sidewall would have the effect of increasing the effective expansion ratio of the 
nozzle and shifting peak performance to a higher NPR, while the twin engine configuration 
would cause interaction effects between the two nozzle jets. Because of the complicated jet 
plume structure at above-design underexpanded NPR's, interaction effects are intense in 
this range and can result in a dynamic "coupling" of the two jets, in which case the twin jet 
configuration behaves much differently than two isolated single jets. Though this 
mechanism generally increases thrust efficiency, it can result in stability and aeroacoustic 
fatigue problems. 
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Conclusions 

In this paper, the basic flow physics and thrust performance characteristics of 
convergent-divergent exhaust nozzles were discussed. From this point, thermodynamic 
relations governing nozzle internal performance were derived using one-dimensional gas 
dynamics, and an approximate integral momentum boundary layer method was developed 
based on the classic Karmdn-Polhausen solution to calculate boundary layer development 
and skin friction losses. Together with angularity effects, these models were written into a 
computational code (NPAC), and used to predict the static thrust performance of 2D-CD 
nozzles. To demonstrate the prediction method, computational results were compared to 
experimental data in four case studies. Important conclusions are as follows. 

1 . The NPAC method showed excellent agreement with experimental data in predicting 
the peak gross thrust efficiency of 2D-CD nozzles for basic 2D geometries (three out 
of the four cases studied to within 0. 1 percent). This leads to an important follow-on 
conclusion relating the physics of nozzle internal flows to the physics represented in 
the NPAC model; namely, losses in 2D-CD nozzles are due primarily to skin friction 
and angularity. Two-dimensional nozzles are generally acknowledged to have lower 
p)eak gross thrust efficiencies than axisymmetric nozzles of the same expansion ratio 
and throat area, but in the past, the general hypothesis has been to attribute these 
additional losses to the existence of comer flows. The current analysis suggests that 
these secondary flows, if they exist, have a minor effect on thrust performance, and 
indicates that lower peak gross thrust efficiencies seen in 2D nozzles are most likely 
due to an increase in wetted area (for the same ID flow area, a square duct has 
approximately 1.13 times more wetted perimeter than a circular duct). 

2. For low expansion ratio nozzles (A9/Ag=1.09, 1.25) with shallow divergence angles 
(a=1.2°, 5.4°), the NPAC thrust performance prediction showed excellent agreement 
with experimental data throughout the entire NPR range. Differences between the 
predicted gross thrust efficiency and experimental data were generally well within the 
approximate 0.5 percent precision limit of the experimental measurement system. 

3 . The NPAC gross thrust efficiency prediction matched experimental data almost exactly 
for a high expansion ratio nozzle (A9/Ag=1.80) with a steep divergence angle (a=l 1°) 
from underexpanded conditions down to overexpanded NPR's at which internal shock 
induced boundary layer separation began to affect nozzle performance. The separation 
had a notable impact on nozzle pierformance, and resulted in a significant increase in 
off-design gross thrust efficiency over what the thermodynamic model predicted. 
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4 . For a twin jet configuration with cutback outside nozzle sidewalls and an extended 
nozzle inter-fairing, the flow physics and nozzle geometry assumptions represented in 
the NPAC model were not suitable to correctly predict nozzle thrust performance. The 
quasi-2D nozzle geometry and dynamic coupling of the two Jets resulted in different 
thrust performance characteristics than would be predicted for an isolated 2D nozzle. 

From these conclusions, some important points can be made. First, the NPAC model 
can accurately predict the static internal thrust performance of basic, isolated 2D-CD 
nozzles. For cases in which shock induced boundary layer separation dominates off- 
design performance, the NPAC prediction cannot account for separation effects, but is 
valid for externally overexpanded, design, and underexpanded NPR's. Though these two 
points define the basic capabilities and limitations of the NPAC method, they in no way 
impose limitations on its application. As was seen in this discussion, one of the most 
powerful uses of the NPAC method lies in interpreting experimental data. By fully 
understanding the fluid mechanics, thermodynamics, loss effects, and physics represented 
in the NPAC model, a comparison of the NPAC prediction with experimental data can lend 
valuable insight into the fundamentals of nozzle thrust performance, and can help explain 
why a particular nozzle might perform as it does. So, disagreement between the NPAC 
prediction and experimental data in some cases may be as valuable as agreement in other 
cases. 

Though NPAC is not meant to compete with more powerful CFD codes, the NPAC 
method is a feasible alternative for nozzle performance modeling when only thrust 
efficiency or boundary layer development predictions are necessary. In some cases, the 
simplicity, ease of implementation, and quick convergence of NPAC may make this a more 
suitable analysis method. In any event, NPAC is a good companion analysis tool to CFD 
codes both as a fundamental prediction method for interpreting CFD results, and also as a 
stand alone, independent fluid dynamic model. 
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